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Abstract 



On ■ 

o\ ■ 
c : 

t-s . We use the closed time-path formalism to calculate fluctuations at phase transitions, 

both in and out of equilibrium. Specifically, we consider the creation of vortices by 
fluctuations, of relevance to the early universe and to 4 ife superfluidity. 

> : 
r- . 

^ : 1 Introduction 

It has become increasingly important to understand the nature of the field fluctuations 
that are a consequence of the phase transitions that occurred in the very early universe. 
On the one hand, fluctuations at the GUT era are assumed to have seeded the large-scale 

Oh! structure visible today. On the other hand, fluctuations at the slightly later electroweak 

o ■ transition are most likely the basis of baryogenesis. 

In this talk I shall consider three aspects of fluctuations. The first is combinatorical, 
establishing a path-integral in terms of which fluctuations can be evaluated. The second 
.£h ! is concerned with the initial- value problem, exemplified by both equilibrium and non- 
equilibrium processes for a real scalar field. The third provides a non-trivial application, 
the creation by fluctuations of (global) vortices in a symmetry-broken U(l) theory of a 
complex scalar field. 

This application is highly relevant, on two counts. Firstly, while there is no unam- 
biguous mechanism for large-scale structure formation in the universe, arguably one of 
the least artificial (originally proposed by Kibble is to attribute it to the spontaneous 
creation of cosmic strings (vortices) by fluctuations at the GUT transitions. [Primordial 
magnetic fields are created by field fluctuations in a related way. See Enqvist, these 
proceedings and elsewhere 0]. 

Secondly, very recent experiments || show that vortices are naturally generated in 
superfluid 4 He on quenching it through its critical density. In the Ginzburg-Landau 
theory of superfluidity the vortices of 4 He are the counterparts of global U(l) cosmic 
strings (albeit non-relativistic). These results have been argued [4| as providing direct 
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support for the creation of cosmic strings in the early universe. We have some observations 
on this. 

This highly enjoyable meeting arrived a little too early for some of the work presented 
here to be completed and my conclusions are rather more qualitative than I would have 
liked. [] A fuller discussion will be given elsewhere, as indicated in the text. 



2 Combinatorics of Fluctuations. 

To be concrete, we begin with the simplest possible theory, that of a real scalar field 0, 
with order parameter ((f)). A natural measure of the fluctuations of 0(i, x) is given by the 
field averages 

(j)(t) = - [ dx(f)(t,x) (2.1) 

V Jv 

over fixed volumes v = 0(l 3 ). [The symbol v denotes both the position of the volume 
in question and its size]. The system as a whole is taken to have a large volume V = 
0(L 3 ) 3> v, in which V is ultimately taken infinite. Since, at any time, the field is 
correlated over some length £, it is sufficient to take I > £. 

To simplify the notation, assume that the possible configurations <&j{x) of in the 
volume V are denumerable (e.g. we impose periodic boundary conditions). Particular 
values of 0„ then single out particular 

In describing the quantum theory of 0, it is convenient to adopt the field representa- 
tion. Suppose, at some initial time t = ti n , the system, described by the state functional 
|$fc,i in ), is an eigenstate of 0, eigenvalue That is, 

0|$ fe ,t in > = $ fe (f)|$ fc ,t iri >. (2.2) 

At a later time t > U n , the state has evolved to 

=$>*il $ i>*> ( 2 - 3 ) 

3 

a superposition of field eigenstates. The probability that the measurement of in | 
gives $j is |cfcj| 2 , most conveniently written in terms of path integrals as 



r<h.{t,x)=<S>j(x) /•<£»(*, x)=*j(x) 

\c kj \ 2 = DfaexpiiStei]} D$ 2 exp{-iS[<h]}. (2-4) 

J <f>l (tin ,x)=<S> k (x) J4> 2 (*in,x)=*fc (x) 



The first term in (2.4) is c k j, the second cjy, where S[(f>] = J d A xC(4>) denotes the classical 
action that determines the evolution of the system. 

The probability p(4> v (t) = 0) that the coarse-grained field 4> v (t) takes some chosen 
value is obtained by checking which $j(x) have their average as 0, and multiplying by 
the probability for finding such a $j. That is, 

P (Mt) - $) = E M 2J (- I <F2H#) - 1) ( 2 - 5 ) 



E / ^01^02 51- d 3 x'Mt,x') - exp{^[0!] - ^[0 2 ]} 

<f>l(tfin)=4>2{tfin) f\ I - „ -\ 

DfaDfa 6 1- d 3 x" 0i(t, x 1 ) - exp{i5[0!] - iS[<j> 2 }}. (2.6) 



<i>\{U n )=<h{Un) = ^k \ v Jv 



1 I have taken advantage of the few weeks since the meeting to develop some of the points a little 
further. Conclusions are unaltered. 
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It might seem that the upper bound of the field integrals in (2.6) should be <pi(t,x) = 
(f>2(t,x), but we can extend the time contour to any final time value tf in > t without 
changing the integral. We shall take tf in infinite. . 

The expression (2.6) is written more simply as the closed-time path integral of Schwinger 
and Keldysh 

p(Mt)~$) = J c D<j)5(^^d 3 x' ( j )l (t,x')-^jexp{tS ctp } (2.7) 

where C = C\ © C 2 is the closed-time path (see below) in which C\ runs from t in to tf in 
and C 2 runs backward from tf in to t in , infinitesimally beneath C\. 



3m(t) 



Co 



Ci 



■ J fin 



- -i(3 



The relevant action S ctp [(f)}, in the presence of a source j(t,x), is described in terms of 
£(0) by 



Sctp[<f>,j] = J c dtdx{C{<P)+j{t)<P{t)) 

/oo 
dtdx (£(00 - £(0 2 ) +h{t)<t> 1 {t) + j 2 (t)<h(t)) 
-oo 

Spatial labels have been suppressed. The doublets are defined to be 

' 4>{t) ted if a = 1. 



*.(*) = 



and 



3a(t) 



[ (f)(t -ie) t - it G C 2 if a = 2. 
j(t) *GCi ifa=l. 



(2.8) 
(2.9) 

(2.10) 
(2.11) 



—j(t — ie) t — ie G C 2 if a = 2. 

Note that, notionally, we have the second leg C 2 of our curve running back e below the 
first. 

This takes account of the quantum uncertainty. However, in general, the initial field 
&(ti n , x) is only specified statistically. Let Pj n [$fc] be the probability that, at time t = ti n , 
the field takes configuration $fc(x). Then 

p(Mt) -$) = T, P ^k] I DfaDfa 5 (- ( d'x'Mt, x') - $) exp{tS ctp } 

k Hl=<t>2=^k \v Jv J 

= [ D<S>P in m I DfaDfo 5 (- ( d 3 x> 0!(t,f) - $) exp{iS ctp }. (2.12) 

J J0 1 =0 2 =ll> \V Jv / 

where we have now relaxed the (artificial) denumerability of the initial state. The gen- 
erality of (2.12) is usually too much for us. To be tractable, simple initial conditions are 
required. The most convenient assumption is that Pj n [$] is Boltzmann distributed as 



P in [$]=Nexp{-p in H in [$]} 
3 



(2.13) 



at some initial temperature Tj n = This permits || us to introduce an initial La- 

grangian density £j n (0) in terms of which Pj n can be expressed as the Euclidean-time 
path integral (the diagonal density-matrix element) 

P m m= I D0 3 exp{-SU0 3 ]} (2-14) 

where Si n [4>] = J dx£i n (4>(x)) is independent of T in . The time contour for the integral 

(2.14) is Cz (see above), running from time t in to U n — i/3i n , over which 3 is periodic. 
We stress that £i n (4>) exists only to parametrise the initial conditions and, in principle, 
has nothing to do with the £{<$) of later times. However, in practice we have in mind a 
situation in which £ in and £ of the closed-time path have the same form, but in which 
the parameters have a time-dependence ||. The other obvious initial condition, in which 
the field <fi is localised at a particular value, in effect, can also be taken as a particular 
choice of Boltzmann distribution, as we shall see. [In principle we can invert (2.14) to 
deduce S in from P in , whether P in describes a Boltzmann distribution or not. In practice, 
it is all but impossible.] 

With this proviso, p of (2.12) can be written as 

p{(j) v {t) - 0) = / DfaDfcDfaS (- I d 3 x'Mt, x) - $) ex P {iS c } (2.15) 

where C is the contour C\ © Ci © C3 and Sc = Si n on C3. 

Let I(x) be the window function (indicator function) for the volume v (i.e. I(x) = 
l,x G v;I(x) = 0,x $ v). On using an exponential representation of the 5-function 

(2.15) , rewritten as 



p(<p v (t) = <p) = J da J DfaDfcDfoexpiiScty] - J dxa(<p - <fa.(t,x))} (2.16) 
can be further recast as 

p((f) v (t) = 0) = J da exp{-iav<p}Z[aI, 0,0] (2.17) 

Z\ji,32,33], defined by 

Z\ji,h,j 3 ] = J DfaDfaDfaexpiiSctyaJa]}. (2.18) 

is the generating functional for Green functions in the presence of sources on all three 
contours (the extension of Sct p [<f) a ,j a ] of (2.8) to include C 3 ). In (2.17) j\ = al(x) 
is a source coupled to the field on C\. The boundary conditions are now 0i(tj n ) = 

h(Un - iflin), 4>l(tfin) = 02 (tfin) 

This splitting into triplets is merely a matter of formal definition. The reason for 
doing this is that it manages to encode the initial conditions into an expression (2.18) 
that looks familiar, and is hence amenable to our usual tricks. For the circumstances 
that we shall consider here, a Gaussian approximation to the fluctuations is a useful first 
step. Adopting the notation 

(F[<f>]) = f DfaDfoDfo F^expiiScfoj = 0]} (2.19) 
** c 

for all F[(j)\, it follows from (2.17) that 

p((f) v (t) = 0) = / daexp{— iav(p}{exp{iav(f) v }). (2.20) 



The Gaussian approximation is obtained by curtailing the cumulant expansion 

(exp{iav<j> v }) = exp{iav(<f) v ) - §aV[(0„0„) - (<f) v } 2 } + 0(a 3 )} (2.21) 

at 0(a 2 ). Performing the a integration gives 

p(Mt) = 0) = iVexp{-072(0^>} (2.22) 

where, for simplicity, we have assumed that (<f> v ) = 0. Not surprisingly, the variance in 
is given by the coarse-grained two-point correlation function. 

More generally, we might wish to calculate the probability p(f v [<f>] — /). where f v [4>] 
is some other coarse-grained functional of 0. In the same approximation 

VUM = /) = Nexp{-p/2(f v f v )} (2.23) 
again assuming (f v ) = for simplicity. 



3 The Initial Value Problem 

As we have observed, the hard work has all been neatly tucked away in the boundary 
conditions. The familiarity of the path integral formalism makes it easy to lose the 
distinction between the contour C3 (the initial condition) and C\(Cq) that carries the 
dynamics. However, for our first case, thermal equilibrium, the distinction has been 
intentionally removed. 



3.1 Thermal Equilibrium. 

This is the most extreme case, in which there is no dynamical evolution of the system. 
Although an unlikely occurrence in the very early universe, thermal equilibrium has been 
examined in great detail, in particular through the thermal effective potential V{4>) 0. 
It is the behaviour of V(<j>) that determines the order of the phase transitions in which we 
are interested. Our approach, given in detail elsewhere || , is close to that of Jona-Lasinio 
||, who defined the effective potential in terms of probabilities for coarse-grained fluctua- 
tions of the type proposed in the previous section. With the probability p(<p v (t) = <p) now 
time-independent, the contours C\ and C2 can be shrunk to nothing, leaving the familiar 
Euclidean-time formalism on contour C3. We stress that, since thermal equilibrium corre- 
sponds to making the choice P» n [$] of (2.14), where Si n [<f>] is now the classical action S[4>] 
of the theory, we have no further freedom in our initial conditions, To cite one common 
misconception, we cannot now use the loop-expanded effective potential following from 
(2.14) as a 'classical' potential in which we subsequently consider the evolution of field 
configurations centred on metastable minima. Such an action corresponds to imposing 
two incompatible boundary conditions simultaneously. However, something different, but 
not wholly dissimilar, may make sense. [See Weinberg flQ| , these proceedings]. 



Let us return to P[<&\ of (2.14), corresponding to equilibrium at temperature T = (3 . 
(With no evolution, all suffices have been dropped). Because of the reality of the integrand 
we can improve upon the Gaussian approximation. The integrated probability satisfies 
the Chebycheff inequality 

P(0„>0) = j D$&v{-PH[$]}0(-J d£ <&{x)-pj (3.1) 

< J D$exp{-f3H[<5>}} exp{(3a (~ J^dx -<^} (3.2) 



for all a > 0. Define uj v (a) by 

exp{(3vu v (a)} = J D<$> exp{-/5iJ[$] + (3 J j$} (3.3) 

where j(x) = al(x), as before. Then 

p(<f>v < 4>) < exp{(3v(u v (a) - a<f))} (3.4) 
an inequality minimised for a = a, 

du„(a) — 



da 

That is 



(3.5) 



p{<t>v < 0) < exp{-(3vVS)} (3-6) 

where V v (4>), the Legendre transform of u> v (a), is the coarse-grained effective potential. 

More specifically, when « m we identify H[§] with ^[^j, the 3-dimensional action 
for the light Matsubara mode, obtained by integrating over all heavy modes (see Kajantie 
TT[ , these proceedings). Thus, but for light-mode self-interactions, V v ((p) gives the usual 



thermal effective potential V((f)) (the convex hull of the loop-expanded potential) in the 
large-f limit, when the inequality (3.6) becomes an equality. We note that, if we were 
considering a gauge theory, rather than just this simple scalar theory, the interesting 
features of the effective potential (e.g.the order of the transition) are largely determined 
by the light gauge field modes, and hence are correctly incorporated in V v . 

We shall be much simpler, initially restricting ourselves to a free scalar field, mass m, 
temperature T ^> m. If I(k) is the Fourier transform of I(x) (normalised to 1(0) = v = 
0(l 3 )) then 

(<p v <Pv) = ^ J(fk\I(k)\ 2 G(k) (3.7) 

where 

G(k) j^— (3.8) 
k z + m z 

is the free-field thermal propagator in this regime. 

What we see from this is that, as we might have anticipated, the effect of coarse- 
graining is, through |J(A;)| 2 , to impose a cut-off at \k\ < ((p v ,(p v ) decreases from 
T/m 2 v at large v (I(k) -^8{k)) to T/Am 2 v = mT/A with A ~ 10, when I = £ = 0{m~ r ). 
That is, the variance in is 

(A0) 2 ~ mT/A (3.9) 

for correlation-volume fluctuations. The value of A depends a little on the shape of V. 
What matters is that it reduces (A0) 2 by an order of magnitude from our first guess |12 



This provides a very useful guide as to when fluctuations are important. [We note that 
V v (4>) gets steeper as v diminishes]. 

We are interested in phase transitions, and for a real scalar theory the only possibility 
is in the breaking of reflection invariance <ft ~ * ~<P- Let us now consider the case when 
the scalar field theory has broken symmetry when cold. We approximate H[$>] by the 
one-loop form 

H[$] = Jdx [|(V$) 2 - > 2 (T)f 2 + ^A$ 4 ]. (3.10) 
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The correlation length £ = m(T) 1 , where m(T) is the effective mass 

m 2 (T)=m 2 (l-?pj (3.11) 

diverges at the critical temperature T c , at which the theory undergoes a second-order 
transition. 

It is not possible to calculate the coarse-grained V v (<f>) of (2.29) analytically. In the 
symmetric phase (T > T c ) perturbation theory is appropriate and we can expand about 
the Gaussian to include a quartic term A„0 4 in the coarse-grained potential V v (4>). Its 
strength X v is, not surprisingly, given in terms of the coarse-grained four-point correlation 
function (<j) v <p v (j) v (j) v ) . Very near to the phase transition the Gaussian approximation 
breaks down, but even then it remains a useful qualitative guide. 

In the symmetry-broken phase perturbation theory breaks down, and even the Gaus- 
sian approximation is difficult. A Gaussian bound can be made by weakening the inequal- 
ity (3.6) further by restricting the spatial integral to x e v in (3.10). Without worrying 
about the details the end result is the highly sensible bound || 

(A0) 2 = (0^> < a 2 + — L- (3.12) 

pvm A l 

The first term on the right hand side of (3.12) describes the long-range field corre- 
lations in the symmetry-broken phase, the second the fluctuations about the minima at 
|0| = a. For correlation volumes v = 0(£ 3 ) = 0(m(T)~ 3 ) the free field result suggests 
that the inequality is improved by replacing Tm(T) = (f3vm 2 (T))~ l by Tm{T) / A. When 
Tm(T)/A = 0(a 2 (T)) there is a significant probability that correlation-volumes of false 
vacuum can be created. This is more transparent if written as the constraint on the 
effective thermal coupling strength A3, that 

A 3 = XT/m(T) = 0(1) (3.13) 



This condition is seen to be just the Ginzburg criterion |L3|], the condition that the 
effective thermal coupling A3 is becoming large, and signalling that the one-loop form 
(3.10) is beginning to break down. [The failure of the Gaussian approximation in the 
symmetric phase is now seen to correspond to | A3 1 =0(1)]. 

If H[§] describes a theory experiencing a first-order transition, as the EW transition is 
expected to do, the situation is more complicated. There has been an extensive discussion 
of the relative role of fluctuations versus quantum tunnelling by Gleiser, Kolb and others 
|Tj| (although our earlier caveat about incompatible boundary conditions is valid here). 
The consensus that thermal fluctuations are inadequate to populate the symmetry-broken 



vacuum state is largely confirmed by our analysis, which will be published elsewhere [15 . 



3.2 Non-Equilibrium Behaviour. 

The picture that we have given above, in which correlation volumes of false vacuum are 
produced with significant probability in the Ginzburg regime demonstrates how large 
fluctuations can occur. However, it may be positively misleading. The expansion of 
the very early universe is so rapid that particles are separated from one another before 
they have undergone sufficient collisions to attain it. To exemplify the nature of out-of 
-equilibrium behaviour, and to offer the strongest contrast to thermal equilibruim, we 
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consider the case of rapid quenching of the scalar theory from an initially symmetric 
state at time t = ti n to a symmetry-broken state, described by the classical action 

S[<f>] = J d A x [±(<90) 2 + > 2 2 + ^A0 4 ]. (3.14) 

More details are given by Lee (these proceedings, and work by Lee, Boyanovsky, de 
Vega |TB[ ) in which a similar calculation is performed, and from which we borrow re- 
sults. [Although motivated by early universe considerations, for simplicity we stay in flat 



spacetime. For a more realistic temporal evolution in the early universe, see de Vega flTTfl 
(these proceedings).] 

At t — ti n we would like the initial state to be symmetric, with Pj n strongly peaked 
about = 0. Despite our caveats about an inability to define temperature in strongly 
non-equilibrium environments our earlier observations suggest that this is most simply 
achieved by taking a Boltzmann probability distribution Pj n [$] = Aexp{ — f3 in H in {Q]}. 
H in [Q] is derived, as in (2.14), from an action 

SM = fdx [\{d<pf - im 2 n 2 ] (3.15) 



for some rrii n . While Si n is a phenomenological action that induces the required initial 
peaking, the effect is the same as if the system were in thermal equilibrium at temperature 
T = T in = /3~^ for all time t < t in . The dispersion in at this initial time is given by (31) 
as (A0) 2 = 0(m in T in ) on coarse-graining, chosen as we wish. The parameter (3 in serves 
solely to fix this. For t > t in the evolution of the initial state is taken to be determined 
by S[(j)] of (3.14 ). We take t in = for convenience. 

In effect, we are taking the theory to be determined for all time by the action 

S t [<f>] = J d A x [\{d<pf + \m\i)<f + ^A(t)0 4 ] (3.16) 

where 

m \t) = { m ™ 2 l <0 n (3.17) 

and 

*(«> = { St ^JJ p is) 

subject to the condition that, for t < 0, it was in thermal equilibrium. In this extreme case 
it is apparent that the equilibrium thermal effective potential V and its coarse-grained 
descendents V v have no role to play and it makes little sense to talk about the order of 
the transition. 

We could have taken A constant for all time, but have not done so in order to delineate 
between its role in establishing the initial condition (for which it is an artefact) and its 
role in the subsequent dynamics, in which it is ultimately a constant of the universe. 
Unsurprisingly, taking interactions into account properly is hard. In practice, because 
of the failure of perturbation theory, it is difficult to do better than a self-consistent 
(Hartree) linearisation of (3.16). Despite our preamble, for our purposes it is sufficient to 
set A(t) = in (3.18) (whence the classical potential is an upturned parabola for t > 0) 
but only to consider the evolution of the system to times at which the field configuration 
would have spread from <fi = to the point of inflexion <fi in f = 0(m/y/\). Since the 
characteristic growth rate is <f>/m = e mt , this means times for which t = 0(m~ 1 /n(l/A)). 
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The spinodal dispersion in <p is now given by the 'free-field' equal time correlation 
function (explicitly dependent on t) 

(<f> v <f> v ) = -j f dxdx 1 (0(i,f)0(t,f)) (3.19) 

= — [ dxdx' I(x)I(x)G(x-x;t) (3.20) 
v l J 

= 4 / d% \I(k)\ 2 G(k;t) (3.21) 
tr J 

in terms of the Fourier transform 

G{k; t) = [ dxG{x; t)e~ i%£ . (3.22) 



At equal-time all correlation functions are the same, once the boundary condition G(x; ti n ) 
G(x; ti n — iftin) nas been implemented. 

G(k; t) possesses both oscillatory modes and exponentially growing (fading) modes, 
as solutions to 

{d 2 + k 2 + m 2 {t))U k {t) = (3.23) 
from which it is constructed. The dominant contributions to (3.21) are the exponential 



modes, and if we take them only then [16 



G(k;t) ~ 9{m * } (l + | ( m 2 +m J; ) (cosh(2ty/m*-fr) - 1) ) coth[\(5 in uj in {k)] 

(3.24) 

where oOi n {k) = \Jk 2 + m 2 n . If we consider fluctuations coarse-grained to volumes v = 
0(l 3 ) the momentum is cut off at \k\ = 0{l~ l ), whereas the ^-function imposes a cut- 
off at \k\ = m. Since we are only interested in I > mT 1 the ^-function is irrelevant. On 
converting to spherical polars in (3.22) k 2 G(k : t) has, for large values of mt, a sharp peak 
at k 2 = 0(m/t). For this to be inside the integral, so that the probability is significant, 



requires I 2 < t/m. Thus the field is correlated on a scale £ = 0(Jt/m). At the relevant 
time t = m _1 /n(l/A) at which domain formation rapidly slows to a halt, the domain size 
is then 

£ = 0(mr 1 (ln(l/\))%) (3.25) 

This is in contrast to £ = 0(m~ l / \f\) in the equilibrium case at the Ginzberg temper- 
ature (3.13). In practice, unless A is extremely small the sizes can be comparable when 
prefactors are taken into account properly. Further details are given by Lee. 



4 Vortex Creation from Fluctuations. 

Consider a field theory invariant under a group G. Suppose the effect of spontaneous 
symmetry-breaking is to reduce G to its subgroup H. The vacuum manifold M can be 
identified with the coset space G/H. If M has a non-trivial first homotopy group rii(M) 
(i.e. the embeddings of loops in M), then vortices (strings) can form |l|j. [If the second 
homotopy group U 2 (M) is non-trivial we can have monopoles, and if the third homotopy 
group n 3 (M) is non-trivial we can have textures]. As we said earlier, the production of 
vortices in the early universe has the capacity to seed large-scale structure through their 
gravitational effects. In general, GUT vortices are assumed to arise from the breaking of 
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a local gauge theory. However, because of their relative simplicity, a considerable effort 
has gone into understanding global vortices, and it is they that we shall mainly consider 
here. Given that the vortices of superfluid 4 He are global vortices, they are of more than 
theoretical interest. 

The simplest theory permitting such vortices is the global U(l) theory of a complex 
scalar field 0, with Minkowski space-time action 

S[<f>] = J d 4 x \\m\f + |m 2 |0| 2 - ^A|0| 4 ]. (4.1) 

With m 2 > the U(l) symmetry is spontaneously broken, and the classical vacuum 
determined by |0| 2 = a 2 = m 2 /X. . The vacuum manifold M = S l has homotopy group 
n 1 (S' 1 ) = Z, and vortices can exist with integer winding number n, Vortices with winding 
number |n| > 1 are unstable to decay into vortices with \n\ = 1. 

Consider an open surface S with oriented boundary dS. The line integral 

measures the winding number of the field configuration on S at time t. (i.e. 27rNs is the 
change in phase of the field as it is taken once around dS). 

If Ns 7^ 0, continuity of requires that it vanish at some point or points of S. A vortex 
of winding number n is a tube of 'false' vacuum (0 ~ 0) containing a line of zeros of for 
which Ns of (4.2) has value n for any closed path enclosing this line. For a large loop dS, 
Ns measures the net winding number of the vortices that pass through it. At a distance 
r from the vortex centre |0| approaches its vacuum value a as |0| = <r(l — 0(e~ mr )). In 
practice it is more convenient to evaluate the related quantity 

N s (t) = f dl.l^jt (4.3) 



-i 



2na 2 
2 



f g dg.(&A 90)^ (4.4) 
f dS!.{d P Ad X ) (4.5) 

J S 



2tt<t 

where, in (3.5), we have used the radial/angular decomposition 

= pe ix (4.6) 

of the field. For a large loop dS the difference between Ns and Ns (not integer) is 
vanishingly small if no vortices pass close to dS, and N s remains a good indicator of 
vortex production. 

On further decomposing the radial field p as p = a + h for Higgs field h, S of (4.1) 
becomes 

S[h, x} = fdx [\{dhf + \a 2 {d X f - m 2 h 2 - \ah 3 - ^\h 4 ] (4.7) 
In the Gaussian approximation the probability that Ns(t) takes the value n is, from (2.23) 

p(N s (t) =n) = exp{-n 2 /2(N s (t)N s (t))} (4.8) 
On further defining the Goldstone mode g by g = ax, from (4.5) it follows that 

(Ns(t)N s (t)) = (^) 2 / j g dS'dS"{{dh' A dg')(dh" A dg")) (4.9) 
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All fields are at time t. The primes (doubleprimes) denote fields in the infinitesimal areas 
dS', dS" of S respectively. For economy of notation we have not made ths scalar products 
explicit. Without loss of generality we take S in the 1-2 plane, whence 

(N s (t)N s (t)) = (^) 2 J J 2 dS'dS"(d i h'd i h"d j g'd j g" - d^dfi' d 3 g' d l9 ") (4.10) 

It is convenient to refine our notation further, decomposing space-time as x = (t, x) = 
(t,XL,xr) where xl = (xi,x 2 ) denotes the co-ordinates of S, and xt = x 3 the transverse 
direction to S. Similarly, we separate 4-momentum p as p = (E,p L ,p T ). . 

Let G h {t,x' - x") = (h(t,x')h(t,x")), G g {t,x' - x") = (g(t,x')g(t,x")) be the Higgs 
field and Goldstone mode correlation functions respectively. As a first step we ignore 
correlations between Higgs and Goldstone fields. That is, we retain only the disconnected 
parts of (NsNs). Eqn. (4.10) then simplifies to 

(N s (t)N s (t)) = (-^) 2 / / dS'dS"[(d i h'd i h")(d j g'd j g") - (a i h'd j h H )(d j g'a iS r)] 

lll(J J Ji,j=1.2 

(4.11) 

which can be written as 

(N s (t)N s (t)) = (^) 2 / fifptfp" G',(p^^)G' fl (^^0|/(pl-pl)| 2 [(pl) 2 (pl) 2 -(pl.pl) 2 ] 

(4-12) 

In (4.12) I{pl) is the Fourier transform of the window function I(xl) of the surface S 
(i.e. I(x L ) — 1 if xl £ S, otherwise zero). The G(p,t) are defined as in (3.22). 

We coarse-grain in the transverse and longitudinal directions by imposing a cut-off in 
three-momenta at \pi\ < A = for some /, as before. Thus N s is now understood as 
the average value over a closed set of correlation-volume 'beads' through which dS runs 
like a necklace. However, we leave the notation unchanged. 

For large loops dS, I{c[l) enabling us to write 

(N s (t)N s (t)) = (JjLyJufr Jd 3 pG h (p+q T ,t)G g (p,t) /#?L|/(&)| 2 [(ft) 2 (?L) 2 -(ft.?L) 2 ] 

(4.13) 

By p + qr we mean (pl,Pt + Qt)- The dependence on the contour dS is contained in the 
final integral 

J = J^ L \I(q L )np L )\q L ) 2 -(p L .q L ) 2 ] (4.14) 
^ irpl J A dq L ql\I(q L )\ 2 (4.15) 

If this is evaluated for a circular loop of radius R, we find 

J = p 2 L 0(2nR/l) (4.16) 

as we might have anticipated. The rms winding number behaves with path length as 
An oc J~2 = O(La) , where L is the number of steps of length /. This is consistent with 
the caorse-grained field phases of different volumes v being randomly distributed. 

The final step is to relate the magnitude of the fluctuations in winding number Ns 
to the magnitude of the radial (Higgs) field fluctuations and angular (Goldstone) field 
fluctuations. To see this requires a specific choice of initial conditions, and we repeat 
those of the previous section, thermal equilibrium and instantaneous quenching. 
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4.1 Thermal Equilibrium. 

To estimate the thermal fluctuations in N$ in equilibrium we a) neglect the positiv- 
ity of p and the Jacobean from the non-linear transformation (4.6) and b) the non- 
singlevaluedness of x- While valid for small fluctuations around the global minima this 
can only be approximate for large fluctuations. With this proviso, at temperature T the 
equilibrium propagators are time-independent, read off from (3.7) as 

°h(P) = -2 ■ T 2 m ( 4 - 17 ) 
p 2 + m z H (T) 

GM = ^ (4.18) 

where m#(T) = \f2m{T) is the effective Higgs mass at temperature T. If we take 
•p\ = 2p*/3 in the integral then, up to numerical factors, we can approximate (3.10) as 

(N S N S ) ~ J (-?—) m(T)T I ctpG h (p) (4.19) 

\2na 2 J J\p\<rn(T) 

where we have coarse-grained to the Higgs correlation length £ = mJ^iT). [We have 
further assumed that the qr integration can be approximated by setting to zero in 
the integrand. Qualitatively this is a reasonable simplification]. The integral in (4.19) is 
essentially the integral (3.21). The end result is that, after substitution, 

(NsNs) = O [ ('^^ )' | 0(L) (4.20) 

where L is the length of the path in units of £. Equivalently, on using our previous results 
for equilibrium 

(N S N S ) = O ((^) 2 ) 0(L) (4.21) 

One power of (h v h v ) comes from Gh, the other from residual factors. That is, the 
phase fluctuations are scaled by the Higgs fluctuations. In the Ginzburg regime, when 
correlation- volumes of the Higgs field can fluctuate to the false vacuum with significant 
probability, the fluctuations in field phase on the same distance scale are of order unity 
and, from the O(L) term, are distributed randomly. [It is unclear whether the inclusion 
of connected correlation functions - linking the Higgs and Goldstone fields - would change 
the results qualitatively in this regime, and is under examination. The inclusion of four- 
point correlation functions in the real scalar field calculations had no dramatic effect, 
although the circumstances were somewhat different.] This ability not to just produce 
'beads' of false vacuum, but to string them as vortices, is the Kibble mechanism referred 
to earlier a common starting-point for numerical calculations ]I8 . 



A similar analysis can be performed for 4 He, treated as a quantum-mechanical system 
of non-relativistic bosonic point-particles with pair- wise potentials V(xi — Xj). In thermal 



equilibrium it is well-known [19] that the grand canonical partition function for such a 



system is identical to the partition function 

Z = J D$ exp{-5[$]} (4.22) 
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of a complex non-relativistic quantum field $. The action in Euclidean space-time (with 
Euclidean time < r < (3) for particles of mass m and (effective) chemical potential fi is 

s[$\ = |^$t(x) ^Iv 2 -^^^) 

+ l[d 4 xf dV \<f>(x)\ 2 V{x - x')mx')\ 2 (4.23) 

V(x) = 5(t)V(x) is the time-instantaneous generalisation of the two-body potential in- 
troduced earlier. 

Equation (4.23) is exact. If we adopt a local approximation for the two-body forces 
that enables us to replace the final term in (4.23) by A|$| 4 , and integrate out the r- 
dependent modes, the end result is a three-dimensional Landau-Ginzburg theory with 
effective Hamiltonian 

H[$] = fdx [— |V$| 2 - MT)|$| 2 + ^A|$| 4 ]. (4.24) 
J 2m 

like that of (3.10). At this level of approximation, in which Euclidean time has disap- 
peared, the distinction between relativistic and non-relativistic theories has also disap- 
peared (except for the details of /x(T)). But for a slight change of definition, repeating 
our analysis will give a variance (An) 2 in winding number (along dS) proportional to the 
path length and to the radial field fluctuation. In thermal equilibrium there is, therefore, 
a true identity between the vortex fluctuations of a relativistic scalar field of the type that 
could occur in the early universe and those of superfluid 4 ife. In the Ginzburg regime, 
when a random phase distribution is relevant, the vortex density will be large. 

As a further example of scalar fields in thermal equilibrium we observed earlier that 
if ri2(M) is non-trivial we can have monopoles. In a cosmological context they are an 
embarrassment [Q], since they would most likely dominate the energy of the universe if 
produced at the GUT phase transition and allowed to survive. As a final generalisation 
of our result (4.21) for winding number we note that it is straightforward to repeat our 
analysis for the production of global monoples. Consider the simplest case of a vector 
triplet of scalar fields in a broken global 0(3) theory. The variance in the monopole 
charge Q v in a volume v is, in the same approximation as the above, a product of three 



two-point correlation functions. The end result is [20] that, in thermal equilibrium at 



temperature T, for coupling constant A and effective Higgs mass m(T), 

(QM = 0((XT/m(T)) 3 )0((m(T) 3 v) 2 ^) (4.25) 

where m 3 v is the volume v measured in terms of correlation volumes, correlation length 
£ = miT)" 1 . In the Ginzburg regime this is understood as the field taking random 
directions in field space in different correlation volumes. 



4.2 Out-Of-Equilibrium Behaviour. 

Our discussion above has shown that, in thermal equilibrium close to the phase transition, 
field fluctuations are large enough to produce local topological charge (winding number, 
monopole charge) at high density, as anticipated by Kibble. However, this is of little use 
unless the future evolution of the system is able to freeze this charge in. Let us return 
to relativistic vortices. The random distribution of field phases shown above leads to 
macroscopic vortex fluctuations, crossing the system volume V. Numerical simulations 
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may take an initial configuration of vortices based on random phases and then, switching 
to zero (or low) temperature solve for their classical evolution. This freezing in of defects 
goes beyond anything that we have said so far. As with the real scalar field of the previous 
section it is simplest to adopt a quenched approximation in which the closed time-path 
action for t > is S of (4.1). In the light of our previous comments one of the more 
interesting choices would be to take the initial state as one of thermal equilibrium in 
the Ginzburg regime with its large fluctuations and watch to see if the vortices freeze. 
However, a simpler, and probably more relevant, choice is to move instantaneously from a 
symmetric to a spontaneously-broken theory, rather as in (3.16), but with the differences 
outlined below. 

Again we assume an initial symmetric state, centred at <fi = 0. Our first guess might 
be to take field configurations Boltzmann-distributed as Pj n [$] = N exp{—f3 in H «„[$]}, 
where H in [<&] is derived from the free-field action 

S m [<P] = /^[i|90| 2 -XJ0| 2 ] (4.26) 

the U(l) generalisation of (2.38). However, when (4.22) is decomposed into radial and 
angular fields as 

Sin[p, x]=J d 4 x [\(dpf + \p 2 (d X f - \m\J\ (4.27) 

the coupling between them (the second term) cannot be ignored since (p) = 0. In the 
absence of any compelling reason (initial thermal equilibrium in the symmetric phase is 
hardly likely in the universe, although it would be appropriate for superfluid 4 He) we 
may as well make a choice of initial condition that is solvable. In this spirit, we take 

S in \p, 7(\=J d 4 x [\{dpf + \a\d X f - \m\J\ (4.28) 

where a 2 = m 2 / A determines the vacuum in the subsequent symmetry-broken phase. 

This Si n , which characterises the initial probabilities, is to be contrasted to the action 
which determines the subsequent spinodal evolution of the field. A similar problem arises 
for times t > U n . Our first thought, the quadratic (in 0) action with downturned potential 

S[p, x \ = j d'x [\{dpf + \p 2 {d X f + \m 2 p 2 } (4.29) 

with ma = V2m, involves interactions between the Higgs and Goldstone fields after 
the radial/angular field decomposition. In the spirit of the Gaussian approximation we 
replace it by 

S[p, x\ = J d 4 x \\{dp) 2 + \a 2 {d X ? + \m 2 p 2 ] (4.30) 

in which the decoupling between p and x is enforced by putting p 2 = o 2 in the second 
term of (4.29) The identity of the massless x action in (4.28) and (4.30) guarantees that it 
remains in thermal equilibrium at all times. Although S of (4.30) is artificial, it is realistic 
in that we would not expect the x modes to be dominantly exponential, but oscillatory. 
This is in contrast to the quenched behaviour of the radial Higgs field p as it falls off the 
hill, with its exponentially growing long- wavelength fluctuations. The Goldstone field, 
always being massless, possesses the IR singularity for small p, 

G g ($t)*y (4-31) 
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With this 1/p 2 effectively cancelling the p\ in J7", (Ns(t)Ns(t)) behaves as 



(N s (t)N s (t)) ~ j(-^)mT in f <tpG p (p;t) (4.32) 

where L is the path length in units of I which, now, we take to be £ of (3.25). 

As before, the fluctuations in winding number vary as the square root of the step 
number, a consequence of random phase fluctuations, with scale set by the fluctuations 
of the Higgs field at the cessation of domain growth. However, with p positive there is now 
only one side of the hill to roll down. The non- vanishing of < p v > prevents the Higgs two- 
point correlation function being identical to that of the scalar field in (3.21) and (3.24). 
However, with < p v > itself showing exponential growth the effect may be qualitatively 
the same. Details will be given elsewhere pi] . How these domains then aggregate to 
produce the true vacuum is beyond our calculations, as is the subsequent evolution of 
the vortex network. What is interesting is the similarity between the equilibrium result 
(4.21) and (4.33). Whether the vortices are produced in equilibrium or strongly out of 
equilibrium the Kibble conjecture seems substantially correct. 

Global monopoles can be treated in the same way. If the Goldstone modes are allowed 
to stay in equilibrium the fluctuations in monopole charge retain the equilibrium volume 
dependence of (4.25), scaled by the Higgs fluctuations. This is not necessarily the case for 
A He. The form (4.23) for the 4 He $-field partition function relied on thermal equilibrium 
for the particle system to be identical to a field theory. We do not know to what extent 
a non-equilibrium system of point particles can be put into correspondence with non- 
equilibrium field theory. However, if it were possible to take a real-time continuation of 
(4.23) as a basis for non-equilibrium calculations there would be no difficulty in principle 
in adopting the same approximations, although the difference between non-relativistic 
and relativistic thermal field theories would now have to be addressed. 

4.3 Local Vortices 

Our experience of continuous symmetries in particle physics has been that all such sym- 
metries have been made local by the presence of gauge fields. In that sense the global 
vortices that we have been discussing are unnatural. 

The local gauge extension of the U(l) theory (4.1) has action 

S[</>, A] = J d A x \-\f^ + \d^ - leArf + m 2 |0| 2 - ^A|0| 4 ]. (4.34) 

(changing factors of | for convenience). This still permits relativistic vortices, the sim- 
plest candidates for local cosmic strings. [Whereas the non-relativistic global strings are 
the vortices of superfluids, the non-relativistic counterparts of local strings are the vor- 
tices of superconductors]. For e 2 /A <C 1 we have a Type-II theory and the strings are 
approximately global, and our previous results should apply. 

The gauge-invariant expression (4.2) for winding number is still valid, but it is more 
convenient to define it through the gauge field as the line integral 

Ns(t) = ^~ I dl.A (4.35) 

Z7T JdS 
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As before, field fluctuations will create local winding number. In the Gaussian approxi- 
mation its variance is 

(N s (t)N s (t)) = (± )V dl[j dl] (Mt^)Mt,x")) (4.36) 

\Z7T/ JdS JdS 

where x',x" denote the positions of the line increments on dS. In one sense this is 
significantly simpler than its scalar counterpart (4.10) since we do not have to worry 
about disconnected and connected parts. The difficulty lies in the more complicated 
correlation function. There is no intrinsic problem in calculating (N$Ns) in thermal 
equilibrium. At the same level of approximation as before (Ns(t)Ns(t)) = 0(e 2 LT), 
up to infrared logarithms. This shows the usual O(L) behaviour and, in step lengths 
£ = m" 1 = (ecr(T)) -1 , (Ns(t)Ns(t)) shows large steps in winding number in the Ginzburg 
regime. Further, since the winding number is essentially the magnetic flux a similar 
calculation can be performed for estimating primordial magnetic fields 0. However, as 
of this moment the work is not complete, but we hope to give the results later. This 
seems a good place to stop. 
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